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ABSTRACT 

We apply a recently developed theoretical model of helium emission to obser- 
vations of both the Orion Nebula and a sample of extragalactic H II regions. 
In the Orion analysis, we eliminate some weak and blended lines and com- 
pare theory and observation for our reduced line list. With our best theoretical 
model we find an average difference between theoretical and observed intensities 
(-f predicted /-^observed — 1) = 6.5%. We argue that both the red and blue ends of the 
spectrum may have been inadequately corrected for reddening. For the 22 highest 
quality lines, with 3499 A < A < 6678 A, our best model predicts observations 
to an average of 3.8%. We also perform an analysis of the reported observational 
errors and conclude they have been underestimated. In the extragalactic analy- 
sis, we demonstrate the likelihood of a large systematic error in the reported data 
and discuss possible causes. This systematic error is at least as large as the errors 
associated with nearly all attempts to calculate the primordial helium abundance 
from such observations. Our Orion analysis suggests that the problem does not 
lie in the theoretical models. We demonstrate a correlation between equivalent 
width and apparent helium abundance of lines from extragalactic sources that 
is most likely due to underlying stellar absorption. Finally, we present fits to 
collisionless case-B He I emissivities as well as the relative contributions due to 
collisional excitations out of the metastable 2s 3 S term. 

Subject headings: atomic data — hii regions — ISM: atoms — ISM: clouds — plasmas- 
Orion — extragalactic 



1. Introduction 



In an era of precision cosmology, accurate theoretical calculations of He I emission are 
essential. In order to determine the primordial helium abundance to a relative accuracy of 
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better than one percent, which is needed to place meaningful constraints on Big Bang Nu- 
cleosynthesis models, He I emissivities also should be known to within a percent. See Bridle 
et al. (2003) for a discussion of the need for accurate measurement of cosmological parame- 
ters. A number of authors have discussed the errors involved (Peimbert et al. 2003, Olive & 
Skillman 2004, Izotov & Thuan 2004). In previous papers, we presented improved calcula- 
tions, in the case-B approximation (Baker & Menzel 1938), of He I emissivities for a range 
of temperatures and densities (Porter et al. 2005, hereafter Paper 1) and for the collisionless 
case (Bauman et al. 2005, hereafter Paper 2). In these papers, we predicted emissivities 
that differed significantly from the previous calculation of He I emissivities (Benjamin et 
al. 1999, hereafter B99). Here we examine the consequences of these new emissivities upon 
abundance determinations. 

As a test of the model helium atom, we predict emission and compare our results with 
the Very Large Telescope observations (Esteban et al 2004, hereafter E04) of the Orion 
Nebula. We first use the simple case-B approximation (in a single, constant temperature, 
homogeneous zone) discussed in Papers 1 and 2, with parameters taken from E04. (See 
Appendix |A] of the present work for fits to both the recombination-only emissivities, as 
a function of temperature, and the collisional contributions due to excitations from the 
metastable 2s 3 S term, as a function of electron density.) As expected, there are some large 
differences, but most of these can be readily explained by optical depth effects not included in 
our simple model. We improve the agreement significantly by selecting a small set of lines and 
optimizing our model by varying helium abundance, temperature, and density. Agreement is 
further improved with a complete simulation, in which the case-B and constant temperature 
constraints are removed and the emission-line region is modeled as an extended region, with 
full radiative transfer effects. In this last comparison we find agreement for 22 highest quality 
He I lines to an average of 3.8%. 

Having validated our model helium atom via comparison with high-quality spectra of the 
Orion Nebula, we then turn to applying our model to primordial helium. We investigate the 
observations of blue compact galaxies by Izotov & Thuan (2004). We find large systematic 
differences in the Izotov & Thuan singly ionized helium abundances determined from different 
lines. We suggest that the systematic differences are primarily due to underlying absorption, 
and we introduce a method by which the uncertainties involved in correcting for this effect 
can be minimized. We discuss other systematic effects and conclude with a summary of our 
results. 

The model helium atom used here is a part of the plasma simulation code CLOUDY 
(Ferland et al. 1998) and is described in Paper 1 and Paper 2. For the purposes of this paper, 
we resolve all nLS terms up to n < 40 and include a series of "collapsed" n-resolved levels 
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with 41 < n < 100. These collapsed levels allow for simpler, less CPU-intensive treatment 
of very highly excited levels and are discussed in detail in Paper 2. At the finite densities 
considered here, this approximation introduces negligible uncertainty. 

2. The Orion Nebula 

2.1. Analysis of Observations 

The E04 Very Large Telescope observations of the Orion Nebula include 100 helium 
emission-line identifications. We dispute the identification of A7937 (identified by E04 as the 
helium line 27d 3 D — 3p 3 P) because of the absence in their observation of lines nd 3 D — 3p 3 P 
with 19 < n < 26 and suggest that the correct identification is more likely the Fe I line at 
A7937.13 (NIST Atomic Spectra Database, version 

3.0.20). 

We note that E04 identify four 
other lines as Fe I emission, although three of these are marked with a question mark to 
indicate an uncertain identification. 

In Figure [TJ we plot the E04 reddening-corrected intensities (relative to H/3) of all lines, 
grouped in series and as a function of the principal quantum number of the upper level. 
A smooth, monotonic progression of these intensities is expected. A cursory examination 
of the trends suggests that many of the weaker lines may be uncertain by a factor of 2 
or more. For the purposes of this study, we discard all lines with /(A) < 0.001 x /(H/3), 
leaving 32 lines. We also discard A8362 (6p 3 P — 3s 3 S) because of a possible blend with CI II 
A8361.84. Compared to Esteban et al. (1998) and Baldwin et al. (2000), the intensity of 
A5048 (4s X S — 2p 1 P) is anomalously high. There is also a possible blend with Fe II A5048.19 
and/or Fe I A5048.43, as evidenced by the fact that the redshift of the observed feature, if 
regarded as the unblended He I line at A5048, would be many standard deviations larger 
than the mean redshift of the remaining helium lines. We therefore discard A5048 as well. 
The remaining 30 lines are given in column 1 of Table [TJ (The other columns of Table [JJ are 
discussed below.) 

As further evidence that some of the E04 line intensity ratios may have uncertainties 
larger than reported, we compare relative intensities of lines originating from the same upper 
level to their theoretical values. Where optical depth effects are not important the theoretical 
ratios of these line intensities should be extremely accurate. The ratio is as follows: 

IjnLS - n'L'S) AjnLS - n'L'S) \" 

I{nLS - n"L"S) ~ A(nLS - n"L"S) A 7 ^ 



1 see |http://physics.nist.gov/PhysRe£Data/ASD/index.html 
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where A' and A" are the wavelengths of the transitions nLS — n'L'S and nLS — n"L"S, 
respectively. In Table [2], we report the theoretical and observed ratios of all pairs of lines 
with the same upper level. The percent difference between the theoretical and observed 
ratios is given, and, in the last column, we report error propagated from the errors given 
by E04. If we assume that the errors reported in E04 represent one standard deviation, we 
would expect (for a normal distribution) about one pair out of the 16 pairs of lines (5%) to 
deviate by more than two standard deviations, but we find that three pairs do (19%). Five 
or six pairs (33%) would be expected to deviate by at least one standard deviation, but we 
find that seven pairs do (44%). Ten or eleven pairs (68%) would be expected to deviate by 
less than one standard deviation, but we find that only nine do (56%). Little changes if we 
exclude the three lines (in two pairs) in Table [2] for which E04 do not report an error except 
to say that it is likely over 40%. Note in particular that three of the pairs deviate by more 
than four standard deviations, which violates the Chebychev Inequality applicable to any 
statistical distribution with finite variance (Hamilton 1964). This analysis suggests that it 
is likely that the errors are underestimated in E04. 

In an attempt to confirm that our 30 lines are correctly identified and not significantly 
blended, we plot in Figure [2] heliocentric recession velocities versus wavelength. We exclude 
A3889 (3p 3 P - 2s 3 S) from the plot but not from our line list because of strong blending with 
H8. The error bars correspond to the least significant digit in the observed wavelengths, as 
reported by E04. The two lines with the greatest recession velocities are A5016 (3p 1 P — 2s l S) 
and A4922 (4c? 1 /} — 2p 1 P), which may be inaccurate due to being near the edge of the 
A4750 — 6800 and A3800 — 5000 wavelength intervals listed in Table 1 of E04. The average 
recession velocity is 15.8 km s _1 , and the standard deviation is 1.4 km s _1 . O'Dell (2001) 
reports 17.9 ± 1.3 km s _1 as a typical heliocentric velocity of "medium ionization" ions, 
which include He + . The low end of O 'Dell's reported range overlaps with the high end of 
the range we find here. We conclude that there are no misidentifications or unknown blends 
in our line list. 



2.2. Theoretical Modeling 

We calculate with CLOUDY a single, homogeneous, constant temperature, constant 
density, case-B model (Model I) using the helium abundance, density, and temperature de- 
rived in the E04 best-fit analysis. In columns 2 of Tabled] we present the fractional difference, 
J(predicted)//(observed) — 1, for each line. In column 3 we present the \ 2 value for each 
line. Optical depth effects are important. Among permitted triplet lines, for example, the 
intensities of both A3889 and A3188 (4p 3 P — 2s 3 S) are overpredicted by case-B calculations, 
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suggesting that, in the observed data, both lines are suffering significant self-absorption due 
to the metastability of 2s 3 S (see Robbins 1968). This argument is strengthened by the fact 
that the intensities of A7065 (3s 3 S — 2p 3 P) and A4713 (4s 3 S — 2p 3 P) are underpredicted in 
case B and apparently enhanced at the expense of the previous two lines. (See the Grotrian 
diagrams of the triplet and singlet systems of helium in Figure [31 ) Because the intensity 
of A4121 (5s 3 S — 2p 3 P) is also underpredicted in case B, it appears likely that this trend 
would extend to it and A2945 (5p 3 P — 2s 3 S). However, we cannot confirm whether A2945 is 
suffering self-absorption because it is outside the wavelength range of the observation. 

There is also some deviation from case-B predictions in singlet lines. For example, 
the intensities of A5016 (3p X P — 2s l S) and A3965 (4p *P — 2s 1 S) are both overpredicted by 
case-B calculations. The first and most likely explanation is that, contrary to the case-B 
assumption, the VUV lines to ground A537.0 (3p l P — ls^) and A522.2 (Ap Y P — ls 1 ^) are 
partially escaping the cloud rather than being completely reabsorbed. Emission in these 
lines will necessarily weaken the intensities of AA5016 and 3965 relative to case B. A second 
explanation is that self-absorption effects from 2s X S are important in these singlet lines in 
the same way they are important for their triplet counterparts AA3889 and 3188, as discussed 
above. If self-absorption were important, however, one would expect the observed intensity 
of A7281 (3s x 5 — 2p X P) to be enhanced relative to case B because each absorption of AA5016 
or 3965 will provide a new opportunity to populate 3s 1 S, necessarily enhancing /(A7281); 
we find the opposite. This point deserves emphasis. We calculate case-B intensities of A7281 
that are 25% greater than in the E04 observation, for which they report an uncertainty of 
8%. If A522.2 partially escapes the cloud, /(A7281) will weaken somewhat relative to case B. 
This effect, however, would also decrease emission in A6678 (3d l D — 2p 1 P), which we do not 
see. We believe, therefore, that /(A7281) being overestimated by case B cannot be explained 
by these optical depth effects. We address this issue again below. 

Next, we modify Model I to allow CLOUDY's optimizer to vary the helium abundance, 
density, and temperature in an attempt to minimize the average y 2 between our predicted 
values and the observed values of a small set of carefully chosen lines. We will label the 
new model "Model II". We exclude any lines with upper principal quantum number n u > 5, 
as well as AA5016, 7281, and 3965 because, as discussed above, they are strongly affected 
by AA537.0 and 522.2 escaping the cloud. In order to minimize the uncertainties involved 
in correcting for the optical depth effects in triplet lines, we introduce a novel approach. 
We note that absorption of a photon at A3889 will inevitably lead either to re-emission of a 
A3889 photon or to decay to 3s 3 S followed by emission of a photon at A7065. Either way, 
the total number of photons in AA3889 and 7065 is conserved. Furthermore, absorption of a 
photon at A3 188 will result in one of the following: 
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1. re-emission of a photon at A3 188; or 

2. decay to As 3 S followed by emission of a photon at A4713; or 

3. decay to 4s 3 S followed by emission of a photon at A3889; or 

4. decay to 3s 3 S followed by emission of a photon at A7065; or 

5. decay to 3d 3 D followed by emission of a photon at AA5876, 3889, or 7065. 

Therefore, the sum of photons in AA3188, 3889, 4713, 5876, and 7065 is independent of the 
optical depth T3889 (and T3i 88 as well, since T3889 and T3188 are related by atomic data alone). 
We include the hydrogen line H8 at A3889 to minimize the uncertainty of deblending this 
line from the close helium line. We use this photon sum, referred to as "triplet photon sum" 
in Table [TJ in our optimization and exclude from our optimization each of the individual 
lines included in the photon sum. (Note that in the vast majority of spectra of H II regions, 
A3 188 is not observed). 

We are left with the following seven quantities to be optimized: the intensities relative 
to H/3 of AA6678, 4922, 3614 (5plP - 2s 1 S), 4388 (h^D - 2p 1 P), 4471 (4d 3 D - 2p 3 P), and 
4121; and the triplet photon sum discussed above. We weight each value by an uncertainty. 
For the line intensities, we use the uncertainties given by E04. For the photon sum, we 
calculate the average of the uncertainties of the individual lines, weighted by photon count, 
and obtain 4%. The optimizer finds best-fit values for the helium abundance (y + = He + /H + , 
by number), log electron density (cm -3 ), and temperature equal to 0.0874 ±0.0005, 3. 941^, 
and 9800 ± 900 K, respectively. The values found in E04 are 0.0874 ± 0.0006, 3.95 ± 0.01, 
and 8730 ± 320 K, respectively. Note that our values were determined using only helium line 
intensities relative to H/3 while E04 determined temperature and density from a number of 
diagnostics that did not include helium intensities at all. Figure H] is a contour plot of the 
average \ 2 between the seven predicted and observed quantities as a function of density and 
temperature, with the helium abundance fixed at the above optimal value. The \ 2 values in 
this plot have been scaled so that the minimum is at unity. The temperature is more tightly 
constrained than the density, which has a rather large standard deviation. 

We substitute our new optimal values into Model I to evaluate Model II, the results of 
which are presented in columns 4 and 5 of Table [H The last row of Table [1] contains the 
average \ 2 f° r each model. The average x 2 obtained with Model II is roughly half of the 
average obtained with Model I. 

Finally, we again use CLOUDY, this time calculating a more realistic model (Model 
III). This model is a plane-parallel slab with abundances typical of H II regions (Baldwin et 
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al. 1991; Rubin et al. 1991; Osterbrock et al. 1992) and grain distributions characteristic 
of those found in the Orion Nebula (Baldwin et al. 1991; van Hoof et al. 2004). The 
slab is heated by a star with a characteristic temperature of 39,600 K (Kurucz 1979). We 
set the surface flux of hydrogen-ionizing photons, $(H) = 10 13 cm~ 2 s _1 and include a 
microturbulence velocity field parameter equal to 16 km s -1 , the average recession velocity 
found above. We enforce constant pressure and include the cosmic ray background. The case- 
B constraint is not specified, and radiative transfer effects (including continuum fluorescence, 
line destruction by background opacities, and line optical depth) are treated self-consistently. 
The results of Model III are presented in columns 6 and 7 of Table [1] and discussed below. 



2.3. Discussion 

The average \ 2 achieved with our Model III is dramatically better than with either 
Model I or Model II. However, several problems remain. The four lines in our set with 
the shortest wavelengths are systematically underpredicted by 16 — 17%. One of them, 
A3 188, could possibly be explained by optical depth effects, although this seems unlikely 
since /(A3889) is accurately predicted. The other three lines are AA3448, 3355, and 3297, 
which are in the same series: np *P — 2s S (for n = 6, 7, and 8). For these three lines and the 
rest of the observed np l P — 2s l S series, we plot in Figure [5] the ratio I(n + 1)/I(n) of both 
the predicted and observed intensities as a function of the principal quantum number n of 
the upper level. The ratios of the predicted values are smoothly increasing, as expected, but 
the ratios of the observed values are anomalously high at n = 5, suggesting that the ~ 16% 
difference noted above is due to inaccuracies in the E04 results. 

Blagrave et al. (2006) have performed a thorough study of reddening toward the Orion 
Nebula, using much higher quality, recent observations. They suggest that the linear ex- 
trapolation of the Costero & Peimbert (1970) reddening function that E04 employ may 
significantly overestimate the extinction shortward of A3500. We support this conclusion. 
The new extinction corrections calculated by Blagrave et al. (based on the analytical extinc- 
tion laws of Cardelli et al. 1989) would also work to reduce or eliminate the discrepancies 
between the present results and E04 for wavelengths longward of about A7000 (see Figure 7 
of Blagrave et al. 2006). We suggest, in particular, that the poor agreement we have shown 
for A7281 may be primarily due to inaccurate reddening corrections. The average difference 
predicted /-^observed — 1 of the lines in Table [T] is 6.5%. If we consider only the 22 lines with 
3499 A < A < 6678 A, we find an average difference of 3.8%. 

Our comparison between theoretical and observed helium intensities is limited primarily 
by the uncertainty in the observations. After eliminating some blends and lines with poor 
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SNR, we have demonstrated good agreement between theoretical and observed values of over 
20 helium line intensities in the Orion Nebula. We believe our theoretical model is capable 
of accurately predicting many more line intensities, including lines originating from larger 
principal quantum number. 

3. Primordial Helium Abundance 

3.1. Analysis of Observations 

Recently, Izotov & Thuan (2004, hereafter IT04) presented extensive data from 33 
observations of blue compact galaxies (BCGs), and used a subset of the strongest He I 
lines to estimate the helium abundance of each system. While IT04 also included in their 
final analysis observations from Izotov & Thuan (1998), Izotov et al. (2001), Guseva et al. 
(2003a), and Guseva et al. (2003b), we consider (for simplicity) only the 33 new observations 
of IT04 in the present analysis. 

We examined the values of y + (A4471), y + (A5876), and y + (A6678) reported by IT04 
in their Table 4 and found statistically significant systematic differences. The average and 
standard deviation of each of these quantities are presented in Table [3j The left panel of 
Figure [6] plots these values versus the weighted means, y + (mean), defined in equation 2 
of IT04. The right panel plots the same values versus the simple means. The left panel 
demonstrates how skewed the weighted means are toward y + (A5876) (because of the much 
stronger signal of A5876) but gives the superficial impression that there is much less deviation 
in the y + (A5876) values than in the y + (A6678) and y + (A4471) values. The right panel 
preserves the relative deviations in each set. The weighted and simple means differ by as 
much as 8%. This analysis is similar to an analysis performed by Skillman et al. (1998) on 
the Izotov et al. (1997) dataset. 

3.2. Sources of Uncertainty 

IT04 derived their abundances using the emissivities of B99, and we find trends similar 
to those shown in Figure [6] when using case-B emissivities from CLOUDY. The largest 
differences are between the abundances determined from A4471 and those determined from 
AA5876 and 6678, amounting to about 7 and 5%, respectively. Note that AA5876 and 6678 
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are lines originating from yrasio levels. Theoretical emissivities from yrast levels may be 
systematically underestimated in the low-density limit by ~ 1 — 2% due to the inaccuracy 
of modelling the helium atom with a finite number of levels (see Paper II). This effect is 
negligible at finite densities as low as 100 cm -3 as collisions force highly excited states to 
local thermodynamic equilibrium and is also too small to explain the systematic differences 
found in y + values. Other effects must be involved. 

One possible explanation is the use of inaccurate reddening corrections. IT04 use the 
interstellar reddening function of Whitford (1958). Olive & Skillman (2004) compared the 
extinction laws of Cardelli et al. (1989) and Whitford in their Figure 2, which shows that use 
of the Cardelli et al. extinction law instead of the Whitford extinction law would increase 
the (reddening-corrected) intensity /(A6678) by about 1.5% relative to the /(A5876), thereby 
increasing the helium abundance necessary to produce the A6678 emission (again, relative 
to A5876). This change would nearly eliminate the systematic difference between the helium 
abundances determined from AA6678 and 5876, but the systematic difference between the 
abundances determined from these two lines and A4471 would remain. The uncertainties 
involved in the reddening correction can be minimized by deriving y + values from a pair of 
helium and hydrogen lines with a small wavelength difference. We calculated ?/ + (A6678) by 
referencing A6678 to Ha (A6563) instead of H/3 (A4861). Similarly, we calculated y + (A4471) 
and y + (A5876) by referencing the lines to ITy (A4340) and Ha (A6563), respectively. The 
y + values obtained are less than the corresponding values derived when referencing the lines 
to H/5, with y + (A5876) and y + (A6678) each decreased by 0.6% and y + (A4471) decreased by 
1.4%. These results are not consistent with the changes produced by using the Cardelli et 
al. law instead of the Whitford law, but both changing the reference line and changing the 
extinction law suggest uncertainties of about 1 — 2%. While certainly significant in primordial 
helium calculations, we find no evidence that reddening corrections are responsible for the 
large systematic differences discussed above. 

For a small subset of systems, IT04 dramatically reduce the systematic differences in y + 
values by accounting for underlying stellar absorption, using theoretical absorption equivalent 
widths (EW a ) from Gonzalez Delgado et al. (1999) and making reasonable estimates for any 
EW a not calculated in that work. Their method does not allow for variations in stellar 
population (i.e., they apply the same set of EW a to every galaxy). The need to allow for 
EW a to vary for different targets was demonstrated by Olive & Skillman (2001) and applied 



2 In regular use in physics literature but rarely in astronomy, the terms yrare and yrast (introduced by 
Grover 1967) are comparative and superlative modifications of the Swedish word yr, meaning "dizzy" , and 
refer, respectively, to levels having high angular momentum and the highest angular momentum (for a given 
principal quantum number). 
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to the IT04 data in Olive & Skillman (2004). IT04 find that after changing EW a (\U71) 
from 0.4 A to 0.5 A their minimization procedure results in an 1% increase in the primordial 
helium abundance Y p , yet they report uncertainies of less than 1% on each of their Y p values. 
Allowing variations between galaxies would surely increase this uncertainty further. 

Recently, Fukugita & Kawasaki (2006) reanalysed the entire new set of IT04 observa- 
tions by considering underlying stellar absorption of the helium lines. It is important to note 
that while Fukugita & Kawasaki do significantly reduce \ 2 i n t ne relation between helium 
abundance and metallicity (dY/dZ), they do not actually model the stellar absorption. In- 
stead they introduce a free parameter that serves to modify the equivalent width of a line. 
Contrary to the analysis performed by IT04, Fukugita & Kawasaki use the same absorp- 
tion equivalent width for each helium line, although (again contrary to IT04) they use a 
different absorption equivalent width for each galaxy. Comparison with a theoretical model 
such as that of Gonzalez Delgado et al. would be beneficial for constraining the Fukugita & 
Kawasaki analysis. See Tremonti et al. (2004) for an application of this technique. 

Figure [7] plots the y + (X) values that IT04 determined from AA4471, 5876 and 6678 
(without correcting for underlying stellar absorption of He I lines) versus the equivalent 
width of the line. Lower y + {\) values tend to correlate with low equivalent width. (A 
similar analysis with a different dataset has been performed in Skillman et al. 1998 with the 
same result.) The trend suggests that excluding lines with an equivalent width less than, 
say, 10 A in primordial helium analyses could reduce the uncertainties involved in correcting 
for underlying stellar absorption. Both Olive & Skillman (2004) and IT04 advocate selecting 
targets using an equivalent width cutoff but only with respect to the equivalent width of 
H/?. By applying a similar additional cutoff to the helium lines, one could, in theory, reduce 
the uncertainty due to underlying stellar absorption to an arbitrarily small amount. We 
note that although a similar effect is accomplished by giving less weight to weaker lines, 
as is done by IT04, that method has the disadvantage that it requires sometimes large 
corrections with unknown uncertainties. The simple mean of all 99 y + (A) values in Table 4 
of IT04 is 0.0810 ± 0.0039. If we exclude lines with equivalent width less than 10 A, the 
mean is 0.0826 ± 0.0034, about 2% larger than in the full set. This procedure has allowed 
us to trade large systematic uncertainties for smaller random uncertainties, the latter of 
which can always be reduced again by making more observations. It is important to note, 
however, that, irrespective of underlying stellar absorption considerations, large systematic 
uncertainties still dominate (see, for example, Peimbert et al. 2002). 
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3.3. Discussion 

We believe that the uncertainties involved in calculating the helium abundance in ex- 
tragalactic systems are much larger than are generally stated. IT04 assign less than 1% 
uncertainty to their value of the primordial helium abundance, but the average y + derived 
from A5876 alone is 7% greater than the average derived from A4471 alone. This discrepancy 
is large, statistically significant, and obtained using both B99 emissivities and emissivities 
predicted by CLOUDY. In our analysis of the E04 Orion Nebula observations, we obtain an 
electron density of nearly 9000 cm" 3 . The BCGs observed by IT04 are less dense by one to 
three orders of magnitude. The uncertainties in the theoretical emissivities in such rarified 
environments are dramatically less than in those found in the Orion Nebula. We therefore 
believe that the present Orion analysis suggests that the systematic discrepancy in y + {\) 
values discussed above is not due to problems with theoretical emissivities. We believe that 
underlying stellar absorption is the largest source of uncertainty. 

4. Conclusions 

We have demonstrated agreement between theory and observation (of the Orion Nebula) 
to an average difference of 3.8% for the 22 lines we designate as having the highest 
quality. 

Higher quality observations of the Orion Nebula could improve agreement between theory 
and observations. 

We support the Blagrave et al. (2006) conclusion that extrapolations of the Costero & 
Peimbert (1970) reddening function shortward of A3500 overestimate extinction to- 
ward the Orion Nebula. We believe the corrections are ps 15% too large. We also 
believe, based upon the Blagrave et al. work, that Costero & Peimbert may have 
underestimated extinction longward of about A7000. 

There is a systematic uncertainty in the IT04 y + values derived from extragalactic sources 
that is most likely due to underlying stellar absorption. The uncertainties involved in 
the correction may be minimized by selecting lines with large equivalent width. 

Much higher quality spectra are essential in further attempts to calculate the primordial 
helium abundance to the accuracy required for tests of Big Bang Nucleosynthesis. 

We thank the National Science Foundation for support through award AST 06-07028 
and NASA through NNG05GG04G. We also thank referee E. D. Skillman for his careful 
review of our original manuscript. 



A. Fits to He I Emissivities 



While calculating emissivities via the above equations sidesteps the self-consistency of 
the CLOUDY calculations, analytical expressions for emissivities are useful in efforts to find 
optimal plasma parameters corresponding to particular observations. In Table HI we present 
four-parameter fits to the case-B, collisionless emissivity of 33 of the strongest He I lines. 
The emissivity, 4:Trj\/n e n- iie +, is calculated as follows: 

47rjA = [a + 6(lnT e ) 2 + clnT e + d/hxT e \ x T~ x x 1(T 25 [ergs cm 3 s -1 ]. (Al) 
n e n He + 

The fits in Table H] introduce negligible new errors, accurately reproducing CLOUDY pre- 
dictions to better than 0.03% for 5000 < T(K) < 25000. 

Collisional contributions to these emissivities are calculated in the same manner as 
Kingdon & Ferland (1995). In Tabled we present the parameters for fits to the collisional 
contributions to the emissivities of any line with the given upper level. We note that col- 
lisional excitations directly affect populations of levels, indirectly affect emission, and are 
independent of the lower level of the enhanced line. The collisional contribution, C/R, is 
calculated as follows: 



- = (1 + 3552 tl a55 /n e )- 1 x ^ a % t\ exp( Ci /t 4 ), (A2) 

i 

where t± is T e /10000, and % is an index that varies from 1 to the number of terms used in 
the fit. As in Kingdon & Ferland (1995), terms comprising less than 1% of the total are 
ignored here. To find the total emissivity of a given line, simply multiply the result obtained 
in Equation IA1I by the quantity 1 + C/R obtained in Equation IA2I via parameters listed in 
Table El 

Figures [H] and [9] compare CLOUDY predictions of collisional contributions to the emis- 
sivity of A5876 with the simple fits of contributions due to excitations from 2s 3 S, as a 
function of temperature at a fixed density, and as a function of density at a fixed tempera- 
ture. The figures indicate that the quantities 1 + C/R introduce negligible uncertainties at 
low temperatures and densities, but underestimate the collisional contribution by as much 
as 5% at high temperatures and densities due to collisional excitations from other levels. At 
such extreme conditions the collisional contributions are themselves uncertain by at least 
5%. At conditions typical of H II regions (n e ~ 100 cm -3 and T e m 15000 K), the error in 
the collisional contribution given by Equation IA2I amounts to a fraction of a percent of the 
total emissivity of A5876. 
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Fig. 1. — Intensities of all helium lines observed by E04, relative to H/5, grouped in series, 
and as a function of the principal quantum number of the upper level. 
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Fig. 2. — Recession velocities of our reduced set of E04 helium lines plotted versus wave- 
length. We have not included A3889 because of the strong blend with H8. 
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3 S 3 p 3 D 1 S 1 p 1 D 

Fig. 3. — Grotrian diagrams of the singlet and triplet systems of the helium atom. The 
effective quantum number, n — 5, is as defined in equation 1 of Paper 2. Labels stacked 
near level nl 2S+1 L indicate the wavelengths in angstroms of transitions nl 2S+1 L — n'V 2S+1 L', 
where I' = L' — I + 1. The labels are arranged vertically in the same sense as the n'l' 2S+1 L' 
levels are arranged. 
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Fig. 4. — Average \ 2 between seven predicted and observed quantities as a function of 
hydrogen density, and electron temperature, T e . 
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Fig. 5. — Ratio I(n +1)/ I(n) of intensities of lines in the np *P — 2s X S series, plotted versus 
the principal quantum number n of the upper level, for both the E04 observations and the 
present (Model III) predictions. The observed ratio appears anomalously high at n = 5. 
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Fig. 6. — IT04 values of y + determined from AA4471, 5876, and 6678 versus (a) weighted 
mean y + (mean) from Table 4 of IT04, and (b) a simple mean, denoted by < y + >. 
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Fig. 7. — IT04 values of y + determined from AA4471, 5876, and 6678 versus equivalent width. 
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Fig. 8. — Ratio of collisional to recombination contributions, C/R, to the emissivity of A5876 
at n e = 1000 cm~ 3 as a function of temperature. 
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Fig. 9. — Ratio of collisional to recombination contributions, C/R, to the emissivity of A5876 
at T e = 15000 K function of electron density. 
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Table I. Predicted intensities, relative to FL3, and fractional differences, 
/(predicted) //(observed) — f, and x 2 values between the line intensities calculated in each 
of three CLOUDY models and observed by E04. The last row gives the average x 2 



obtained for each model. 


Wavelength ^AirJ 


Model I 




Model II 




Model III 






I(\)/I(H(i) 


Diff 


V 2 

X 




Diff 


V 2 

X 


iw/i(H0) 


Diff 


V 2 

X 


3188 a 


0.041 


0.456 


32.542 


0.042 


0.520 


42.197 


0.023 


-0.175 


7.055 


3297 a 


0.001 


-0.159 


0.394 


0.001 


-0.146 


0.324 


0.001 


-0.168 


0.454 


3355 a 


0.002 


-0.187 


3.136 


0.002 


-0.176 


2.707 


0.002 


-0.170 


2.466 


3448 a 


0.003 


-0.177 


5.695 


0.003 


-0.164 


4.783 


0.003 


-0.156 


4.239 


3499 


0.001 


-0.051 


0.072 


0.001 


-0.058 


0.095 


0.001 


-0.051 


0.072 


3513 


0.001 


-0.015 


0.008 


0.001 


-0.022 


0.017 


0.001 


-0.015 


0.008 


3530 


0.002 


-0.072 


0.189 


0.002 


-0.078 


0.223 


0.002 


-0.071 


0.180 


3554 


0.002 


0.006 


0.003 


0.002 


0.001 


0.000 


0.002 


0.010 


0.008 


3587 


0.003 


-0.032 


0.135 


0.003 


-0.037 


0.180 


0.003 


-0.027 


0.093 


3614 


0.005 


-0.002 


0.001 


0.005 


0.017 


0.061 


0.005 


0.009 


0.018 


3634 


0.005 


-0.045 


0.458 


0.005 


-0.050 


0.563 


0.005 


-0.038 


0.315 


3705 


0.007 


-0.003 


0.003 


0.007 


-0.008 


0.027 


0.007 


0.008 


0.028 


3820 


0.012 


-0.042 


2.140 


0.012 


-0.048 


2.802 


0.012 


-0.026 


0.816 


3889 b 


0.217 


0.446 


221.455 


0.224 


0.490 


266.437 


0.154 


0.021 


0.513 


3927 


0.001 


-0.006 


0.009 


0.001 


-0.011 


0.027 


0.001 


-0.007 


0.011 


3965 


0.010 


0.055 


3.301 


0.010 


0.077 


6.506 


0.010 


0.023 


0.591 


4009 


0.002 


0.095 


3.590 


0.002 


0.088 


3.119 


0.002 


0.091 


3.329 


4121 


0.002 


-0.158 


22.130 


0.002 


-0.070 


3.507 


0.002 


0.007 


0.033 


4144 


0.003 


0.074 


3.458 


0.003 


0.067 


2.837 


0.003 


0.070 


3.047 


4388 


0.006 


0.037 


3.370 


0.006 


0.032 


2.488 


0.006 


0.032 


2.547 


4471 


0.045 


-0.002 


0.040 


0.045 


-0.002 


0.040 


0.047 


0.033 


10.563 


4713 


0.005 


-0.258 


1206.619 


0.006 


-0.171 


427.810 


0.007 


0.070 


48.490 


4922 


0.012 


-0.007 


0.434 


0.012 


-0.014 


1.990 


0.012 


-0.013 


1.760 


5016 


0.026 


0.121 


146.383 


0.027 


0.150 


223.831 


0.024 


0.040 


15.917 


5876 


0.128 


-0.115 


18.805 


0.129 


-0.103 


14.755 


0.134 


-0.071 


6.575 


6678 


0.035 


-0.093 


2.887 


0.035 


-0.103 


3.697 


0.035 


-0.100 


3.433 


7065 c 


0.034 


-0.544 


289.341 


0.039 


-0.467 


157.203 


0.082 


0.115 


2.694 


7281 c 


0.007 


0.252 


9.904 


0.008 


0.375 


21.997 


0.008 


0.287 


12.879 


10027 c 


0.002 


0.089 


0.311 


0.002 


0.059 


0.135 


0.002 


0.073 


0.209 


10311 c 


0.001 


0.027 


0.028 


0.001 


0.021 


0.017 


0.001 


0.043 


0.072 


triplet photon sum d 


0.409 


-0.042 


1.211 


0.426 


-0.002 


0.003 


0.427 


0.000 


0.000 


average x 2 






63.808 






38.399 






4.142 



a Reddcning corrections may be significantly overestimated. 
b Blended with H8. 

c Rcddening corrections may be significantly underestimated. 
d See text for discussion. 
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Table 2. Theoretical and observed relative line intensity ratios of lines with the same 

upper level. 



n 2S + 1 L 


n' 2S + 1 L' 


n" 2S+l L » 


Wavelengths (Air) 


!'/!" 


DifTcrence 


upper 


lower 


lower 


A' (A) 


A"(A) 


Theor. 


Obs. 


in Std Dev 


6d 3 D 


2p 3 P 


3p 3 P 


3820 


10311 


8.71 


9.3±1.5 


0.4 


7d 3 D 


2p 3 P 


3p 3 P 


3705 


9517 


8.16 


23.9±3.8 


4.2 


7d 1 D 


2p x P 


3p l P 


4009 


10138 


6.68 


6.3±1.1 


-0.3 


7p x P 


2s ^ 


3s x 5 


3355 


8915 


10.10 


11.1±2.2 


0.4 


8d 3 D 


2p 3 P 


3p 3 P 


3634 


9063 


7.86 


9.5±1.5 


1.1 


8p x P 


2s 1 S 


3s J 5 


3297 


8518 


9.69 


13.5±4.8 


0.8 


9d 3 D 


2p 3 P 


3p 3 P 


3587 


8777 


7.67 


4.3±0.7 


-5.1 


llp 3 P 


3s 3 S 


3d 3 D 


7062 


8854 


4.63 


2.4±0.5 


-4.6 


lid 1 !? 


2p 1 P 


3p x P 


3806 


8931 


6.15 


11.0±3.4 


1.4 


12p 3 P 


3s 3 S 


3d 3 D 


6990 


8740 


4.71 


4.3±1.3 


-0.3 


12d 3 D 


2p 3 P 


3p 3 P 


3513 


8342 


7.39 


6. Oil. 3 


-1.1 


12d 1 D 


2p x P 


3p x P 


3785 


8817 


6.10 


7.2±3.3 


0.3 


Ud 3 D 


2p 3 P 


3p 3 P 


3488 


8204 


7.30 


9.7±2.9 


0.8 


15d 3 D 


2p 3 P 


3p 3 P 


3479 


8156 


7.27 


7.8±2.4 


0.2 


l6d 3 D 


2p 3 P 


3p 3 P 


3472 


8116 


7.24 


10.5±3.8 


0.8 


17d 3 D 


2p 3 P 


3p 3 P 


3466 


8084 


7.22 


18.0±10.2 


1.1 
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Table 3. Average and standard deviation of the helium abundances determined from 

Table 4 of IT04. 



A 


<y+W > 


a 


4471 


0.0778 


0.0034 


5876 


0.0835 


0.0033 


6678 


0.0818 


0.0027 
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Table 4. Parameters for fits to He I emissivities, 4nj\/n e n ne +. See text for formula. 



Wavelength (Air) 


n 2S+1 L 


n 2S+1 L 










A 


upper 


lower 


a 


b 


c 


d 


2945 


5p 3 P 


2s 3 S 


-1.1849E+05 


-4.2559E+02 


1.2913E+04 


3.5306E+05 


3188 


4p 3 P 


2s 3 S 


-2.3591E+05 


-8.6438E+02 


2.5969E+04 


6.9653E+05 


3614 


5p 1 P 


2s Y S 


-3.0442E+04 


-1.2890E+02 


3.5599E+03 


8.5434E+04 


3889 


3p 3 P 


2s 3 S 


-6.0645E+05 


-2.5058E+03 


7.0181E+04 


1.7184E+06 


3965 


4p!p 


2s l S 


-6.2095E+04 


-2.6386E+02 


7.2700E+03 


1.7424E+05 


4026 


5d 3 D 


2p 3 P 


-4.8009E+04 


-3.4233E+02 


7.5849E+03 


9.2951E+04 


4121 


5s 3 S 


2p 3 P 


1.0266E+04 


6.5524E+01 


-1.3829E+03 


-2.6154E+04 


4388 


5d l D 


2p 1 P 


-1.1349E+04 


-8.7804E+01 


1.8895E+03 


1.9916E+04 


4438 


5s 1 S 


2p 1 P 


3.7527E+03 


2.2646E+01 


-4.8794E+02 


-9.9429E+03 


4472 


4d 3 D 


2p 3 P 


-3.5209E+04 


-4.5168E+02 


8.5367E+03 


9.1635E+03 


4713 


4s 3 S 


2p 3 P 


2.4264E+04 


1.5802E+02 


-3.2993E+03 


-6.1132E+04 


4922 


4<2 1 D 


2p x P 


-6.1979E+03 


-1.1414E+02 


2.0047E+03 


-8.6195E+03 


5016 


3p x P 


2s 1 S 


-1.3442E+05 


-5.9029E+02 


1.6033E+04 


3.7142E+05 


5048 


4s *S 


2p x P 


7.1813E+03 


4.5389E+01 


-9.5007E+02 


-1.8788E+04 


5876 


3d 3 D 


2p 3 P 


2.0620E+05 


1.7479E+02 


-1.3548E+04 


-7.3492E+05 


6678 


3d x D 


2p 1 P 


6.7315E+04 


7.5157E+01 


-4.7101E+03 


-2.3610E+05 


7065 


3s 3 S 


2p 3 P 


5.8675E+04 


4.1954E+02 


-8.2053E+03 


-1.4565E+05 


7281 


3s 1 S 


2p 1 P 


1.3544E+04 


9.9057E+01 


-1.8831E+03 


-3.4394E+04 


9464 


5p 3 P 


3s 3 S 


-6.5519E+03 


-2.3532E+01 


7.1399E+02 


1.9522E+04 


10830 


2p 3 P 


2s 3 S 


-3.9020E+05 


-1.7846E+03 


4.9448E+04 


1.0443E+06 


11969 


5d 3 D 


3p 3 P 


-4.8419E+03 


-3.4526E+01 


7.6497E+02 


9.3744E+03 


12527 


Ap 3 P 


3s 3 S 


-7.5547E+03 


-2.7681E+01 


8.3162E+02 


2.2306E+04 


12785 


5f 3 F 


3d 3 D 


3.7439E+04 


1.0086E+02 


-3.4590E+03 


-1.2157E+05 


12790 


5f 1 F 


3d x D 


1.2480E+04 


3.3619E+01 


-1.1530E+03 


-4.0523E+04 


12968 


bd^ 


3p x P 


-1.4362E+03 


-1.1111E+01 


2.3910E+02 


2.5205E+03 


15084 


4p 1 P 


3s l S 


-3.3007E+03 


-1.4026E+01 


3.8644E+02 


9.2617E+03 


17002 


4d 3 D 


3p 3 P 


-2.4898E+03 


-3.1940E+01 


6.0368E+02 


6.4809E+02 


18685 


4/ 3 F 


3d 3 D 


8.7712E+04 


2.8178E+02 


-8.8285E+03 


-2.5836E+05 


18697 


4/ X F 


3d x D 


2.9217E+04 


9.3839E+01 


-2.9404E+03 


-8.6063E+04 


19089 


4d x D 


3p x P 


-5.7257E+02 


-1.0544E+01 


1.8518E+02 


-7.9601E+02 


19543 


Ap 3 P 


3d 3 D 


-4.4053E+03 


-1.6141E+01 


4.8493E+02 


1.3007E+04 


20580 


2p x P 


2s 1 S 


-6.5411E+04 


-3.4797E+02 


8.7844E+03 


1.6613E+05 


21118 


4s 3 S 


3p 3 P 


3.7041E+03 


2.4122E+01 


-5.0366E+02 


-9.3323E+03 
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Table 5. Fit parameters for the relative collisional contributions, C/R, to the emissivities 
of lines with the given upper level. See text for formula. 



upper 


i 




bi 




3s 3 S 


1 


37.2702 


-1.2670 


-3.3640 




2 


2.6982 


-1.2918 


-3.6989 




3 


0.9598 


-1.3903 


-4.5122 


3s *S 


1 


17.4945 


-1.4946 


-3.5982 




2 


0.5015 


-1.1724 


-4.5518 


3p 3 P 


1 


8.9027 


-1.0970 


-3.6989 




2 


1.0004 


-1.0491 


-4.3800 




3 


0.5387 


-0.5550 


-4.5449 




4 


0.3533 


-1.1615 


-4.8186 




5 


0.2345 


-0.6632 


-4.9006 


3d 3 D 


1 


6.7937 


-0.1116 


-3.7761 




2 


0.1808 


-0.8306 


-4.5122 




3 


1.3478 


-0.4017 


-4.5459 




4 


0.4792 


-0.4062 


-4.9012 




5 


0.2950 


-0.8224 


-4.9013 


3d^ 


1 


0.4340 


-0.7808 


-3.7766 




2 


0.1942 


-0.7687 


-4.5459 




3 


0.1263 


-1.1108 


-4.9012 




4 


0.0549 


1.8948 


-4.9013 




5 


0.0785 


2.0453 


-5.0942 




6 


0.0935 


2.0461 


-5.0942 


3p x P 


1 


2.8310 


-1.0005 


-3.7917 




2 


0.8019 


-1.3326 


-4.4724 




3 


0.4130 


-1.0703 


-4.5452 




4 


0.2782 


-1.3273 


-4.8641 




5 


0.1912 


-1.0948 


-4.9008 


4s 3 S 


1 


29.1613 


-1.3278 


-4.3800 




2 


1.2121 


-1.4745 


-4.5122 


4s X S 


1 


17.4301 


-1.5851 


-4.4724 




2 


0.3277 


-1.3035 


-4.9043 


4p 3 P 


1 


8.5965 


-1.2591 


-4.5122 




2 


0.6886 


-1.2274 


-4.8186 




3 


0.2282 


-0.7294 


-4.9006 


4d z D 


1 


4.4397 


-0.2954 


-4.5449 




2 


0.1341 


-0.7331 


-4.8839 




3 


0.7546 


-0.5041 


-4.9012 


4d 1 D 


1 


2.5507 


-0.8404 


-4.5452 




2 


0.3106 


-0.8657 


-4.9012 




3 


0.0739 


2.0732 


-5.0942 




4 


0.0399 


2.0304 


-5.0942 


4/ 3 F 


1 


3.1027 


-0.2421 


-4.5459 




2 


0.6790 


-0.6626 


-4.9013 




3 


0.0348 


2.3128 


-5.0942 



-29- 



Table 5 — Continued 



upper 


j 






Ci 




4 


0.0629 


2.3128 


-5.0942 




1 


0.9680 


-0.6405 


-4.5459 




2 


0.2817 


-0.9699 


-4.9013 




3 


0.1045 


2.3128 


-5.0942 




4 


0.1885 


2.3128 


-5.0942 


4p!p 


1 


2.8275 


-0.9962 


-4.5518 




2 


0.5209 


-1.3550 


-4.8641 




3 


0.2236 


-1.1224 


-4.9008 




4 


0.0641 


1.8784 


-5.0730 




5 


0.0385 


1.8788 


-5.0940 


5s 3 S 


1 


33.8477 


-1.4584 


-4.8186 




2 


0.4717 


-1.2851 


-4.8839 


5s *S 


1 


17.8299 


-1.5992 


-4.8641 


5p 3 P 


1 


5.7783 


-1.0686 


-4.8839 


5d 3 D 


1 


4.0162 


-0.4399 


-4.9006 


bd l D 


1 


2.7165 


-0.9041 


-4.9008 




2 


0.0841 


2.0783 


-5.0942 


5/ 3 P 


1 


3.8606 


-0.3120 


-4.9012 




2 


0.0621 


2.2606 


-5.0942 


5/lF 


1 


1.5158 


-0.6943 


-4.9012 




2 


0.1860 


2.2606 


-5.0942 


5p x P 


1 


3.4972 


-1.0841 


-4.9043 




2 


0.1006 


1.8740 


-5.0730 




3 


0.0413 


1.8744 


-5.0940 



